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Abstract 

We show that quantum correlations as quantified by quantum discord can characterize quantum phase 
transitions by exhibiting nontrivial long-range decay as a function of distance in spin systems. This is 
rather different from the behavior of pairwise entanglement, which is typically short-ranged even in critical 
systems. In particular, we find a clear change in the decay rate of quantum discord as the system crosses 
a quantum critical point. We illustrate this phenomenon for first-order, second-order, and infinite-order 
quantum phase transitions, indicating that pairwise quantum discord is an appealing quantum correlation 
function for condensed matter systems. 
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1. Introduction 

A quantum phase transition (QPT) [X] is primar- 
ily characterized by a qualitative sudden change in 
the ground state of an extended quantum system 
as an external parameter or an internal coupling is 
continuously varied. QPTs occur effectively at zero 
temperature and its typical quantum fluctuations 
are believed to be due to genuine quantum correla- 
tions. In recent years, concepts from quantum in- 
formation theory have been extensively applied to 
the study of QPTs in quantum many-body systems 
0, In particular, the existence of quantum cor- 
relations has usually been inferred by the presence 
of entanglement among parts of a system. Indeed, 
entanglement displays a rather interesting behavior 
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at QPTs, being able to indicate a quantum criti- 
cal point (QCP) through nonanalyticities inherited 
from the ground state energy 0, 0]. This behav- 
ior is already observed from pairwise measures of 
entanglement as given, e.g., by concurrence [6| and 
negativity Q- 

Although pairwise entanglement measures usu- 
ally exhibit scaling behavior at a QCP, they are 
typically exponentially short-ranged To some 
extent, long-distance pairwise entanglement may be 
engineered in several many-body systems by conve- 
niently setting microscopic parameters. However, 
this adjustment does not coincide with QCPs in 
those systems nor critical scaling of entanglement 
with distance is observed [Si]. This is somewhat 
surprising since one expects a slow decay of quan- 
tum correlations in a quantum critical region. In 
this Letter, we investigate this problem by quanti- 
fying quantum correlations through quantum dis- 
cord (QD) [9, 10]. In particular, we show that QD 
provides the expected long-range behavior of quan- 
tum correlations for several spin chains exhibiting 
QPTs. Such a behavior is remarkable since, even 
though QD has been shown to be non-vanishing for 
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almost all quantum states llj, its decay pattern 



as a fu notion of distance is unresolved in general 
grounds. We illustrate our results for the transverse 
field XY chain, where long-distance QD is achieved 
due to the onset of magnetic order, as well as for the 
XXZ chain in the presence of domain walls, where 
the long-range behavior of quantum correlations is 
obtained as a consequence of the polynomial decay 
of QD. 

2. Quantum discord 

Before introducing quantum discord, let us call 
up some concepts from classical information the- 
ory (CIT). In CIT the uncertainty about a ran- 
dom variable A, which can assume the values 
a (from a set A) with corresponding probability 
Pa := Pr(74 ~ a), is given by the Shannon's en- 
tropy H{A) = -J2aPa^og2Pa- The Uncertainty 
about two random variables A and B taken to- 
gether reads H{AB) = - Y.a,bPaA^°Z2Pa,b-, with 
{Pa.b '■— Pr(j4 — a,B — b)} being the joint prob- 
ability distribution. The total amount of correla- 
tion between A and B is given by the difference 
in the uncertainty about A before and after B is 
known, i.e., J{A : B) = H{A) - H{A\B), where 
H{A\B) = -Y.a,bPa:b^0E2Pa\b IS the Conditional 
entropy, with Pa\b standing for the probability for 
A — a given that B — b. From Bayes' rule, 
Pa\b — Pa.b/Pb, we Can rewrite J{A : B) in the 
equivalent form I( A : B) = H{A)+H{B)-H{AB). 

In the quantum domain, the uncertainty over 
the state p is given by the von Neumann's entropy 
S{p) — — tr(plog2p). Thus, a straightforward gen- 
eralization of T{A : B) to the quantum realm read- 
ily follows as 



I{Pab) := S{pa) + S{pb) - S{pab)- 



The quantity I{pab) is called quantum mutual in- 
formation and is a well established information- 
theoretic measure of the total (quantum plu s clas- 
sical) correlations in a quantum state [12| . On 
the other hand, the average uncertainty about the 
state of the system A after a complete set of pro- 
jective measurements {11^} is performed on sys- 
tem B is given by Sn{Pa\b) = J2jPjS{pi), where 
Pj = trb(njpf,) and p{ = tn{IijPab^])- Thus, a 
quantum version of J^{A : B), independent of mea- 
surement "direction" in the Hilbert's space, can be 
defined as 



J{Pab) S{pa) 



nYmSn{Pa\b)- 

J 



(2) 



In Ref. [3], OUivier and Zurek noticed that, while 
the classical expressions for mutual information 
given by J{A : B) and I{A : B) are equivalent, 
their quantum generalizations ^ and ([2]) can be 
different depending on the state of the system. This 
difference originated the quantum discord Q{pab), 
which reads 



Q{Pab) := I{Pab) - J{Pab)- 



(3) 



The QD is a non-negative asymmetric quantity 
that vanishes only for states encoding a joint clas- 
sical probability distribution. 



^Pt]Va){ia\ ® \3b){jb\ 



(4) 



where {|ia)} ({Ijb)}) form an orthonormal basis for 
the system A {B) (see Refs. [l3, 14| for a discus- 
sion about the symmetry properties of quantum 
correlation quantifiers). QD is a measure of cor- 
relations based on information-theoretic concepts 
and is intended to capture all quantum correlations 
present in a quantum state [9|. Recently, it has 
received a great deal of attention, exhibiting re- 

, be- 
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markable behavior under decoherence 
ing experimentally tested [2^-[2^ and displaying 
applications in several contexts [27|- 34 [. Concern- 
ing QPTs, QD has been considered as an indica- 
tor of QCPs in several systems, having succeeded 
in this task even in situations where entanglement 
fails [13-0. 



3. The XY model 

In order to investigate the scaling behavior of QD 
as a function of distance, let us begin by consider- 
ing a chain of spin- 1/2 particles anisotropically in- 
teracting in the xy spin plane and subjected to a 
magnetic field in the z-direction. This system is de- 
scribed by the XY model, governed by the following 
normalized Hamiltonian 



H 



X 



1=0 
i=0 



(5) 



with N being the number of spins in the chain, 
(7™ the i-th spin Pauli operator in the direction 
TO — x,y, z and periodic boundary conditions are 
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assumed, i.e., cr]^' — cr™. The anisotropy 7 is con- 
strained to the range < 7 < 1. For 7 — )■ 0, the 
model reduces to the XX model whereas for all the 
interval < 7 < 1 it belongs to the Ising univer- 
sality class, reducing to the transverse field Ising 
model at 7 = 1. The dimensionless parameter A is 
proportional to the reciprocal of the external trans- 
verse magnetic field. For A = 1, a quantum critical 
line takes place for any value of 7 in the range con- 
sidered in this work. 

The exact analytical solution of the XY model in 
the thermodynamical limit {N — (X)) is well known 
(41I . The Hamiltonian ([5]) can be diagonalized 
via a Jordan- Wigner map followed by a Bogoliubov 
transform. By considering the thermal state at fi- 
nite temperature and taking into account the Z2 
symmetry of the XY model, namely, the invariance 
under parity transformation (^j^ i o" ? 1 the reduced 
state for the spins and n reads [43[ 



POr. 



(6) 



where I"" is the identity operator acting on the joint 
state space of the spins and n. Due to the fact that 
the system is translationally invariant, the reduced 
state ^ depends only on the distance n — \i ~ j\ 
between spins i and j. The magnetization density 
(cr^) as well as the two-point correlation functions 
i'^o'^n) '^^^ be directly obtained from the exact so- 
lution of the model [4l|, (their expressions are 
made explicit in Appendix A). We will use through- 
out this Letter the zero temperature limit (T — > 0) 
of the reduced thermal state given in Eq. (|6]), which 
is called thermal ground state. This unbroken state 
was shown to provide an exact description of the 
critical behavior of entanglement as well as its scal- 
ing in finite systems 0,0, Hi- The total informa- 
tion shared by the spins in the state ^ is given 

by 



with 



HfOn) = S{po) + S{pn) ~ S'(pon), 



S{pq) ^ S{pn) 



(7) 



^1±™,,^^1±(^ (8) 



i=0 



and 



S{pQn) ^Y{S.i^Og2ii+r]t\og2m), (9) 



with 

6 = (1 

+(-l)'v/(K<) - WW + 4(a^>^]/4, 
and 

v^ = [1 - (^g^ + {-iy{{^oO + {«))] /4- 

Following j35{ . we numerically verified that the 
minimum in Eq. ^ is attained, for all values 
of A, 7, and n considered in this Letter, by the 
following set of projectors: {|+)(+|, |— )(— |}, with 
|±) = (It) ± |;))/\/2, where {|t) , are the eigen- 
states of (T^. Thus one obtains 



where Hh 



J{pOn) = HbiniPl) + HbiniP2), 

inix) is the binary entropy 



(10) 



Hun{x) = -a;log2a; - (1 - a:)log2(l - x) (11) 



and 



Pi 

P2 



2(1 + ^)). 



(12) 



i=0 



This provides, therefore, an analytical expression 
for evaluating QD for a distant pair described by 
Pon - We note that the algorithm proposed in Ref. 
[45j to evaluate QD in X states is indeed valid in 
our case if Ko-qCt^)! > \{(7Qcr'^)\ 0. Using the an- 
alytical solution for the XY model, one can verify 
this inequality for the whole range of values of A, 
7, and n considered in this Letter. Moreover, as 
shown in [46^, the difference between the QD ob- 
tained using the algorithm of Ref. 45] and that 
obtained by maximizing over all measurements is 
very feeble. Thus, by applying Eas.([7|)- (rT^ . we can 
reliably analyze the decay pattern of correlations in 
the XY chain. 

A change in the range of correlations is a typ- 
ical indication of a QPT in a many-body system. 
Indeed, such a change is shown in Fig. [TJ where 
we plot a pattern of decay of QD as function of A 
and 7. This figure was constructed by computing 
the ratio between the quantities QiPon) and 

MQ{pqi), with M = 10, for each pair (7, A). In 
other words, we compare the actual area under the 
curve Q{pon) with the bare case where QD remains 
constant with distance. Although Fig. [T] itself does 
not reveal the decay rate of QD, it clearly shows a 
qualitative change in the distance behavior of QD 
as the system passes through the critical line A = 1. 
Indeed, a slower decay of QD is exhibited in the re- 
gion A > 1, where magnetic order takes place. This 
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Figure 1: (Color Online.) Pattern of decay of quan- 
tum discord with the distance between the spins 
sites as a function of A and 7. Colors toward the 
red indicate a slower decay of quantum discord with 
distance (see the text for details). 



change in the range of pairwise quantum correla- 
tions at A = 1 is made explicit in Fig. [51 where we 
plotted QD for a spin pair as a function of the dis- 
tance n. The curves in Fig. [5] are the exponential 
fits of quantum discord. We observe that, for both 
examples of anisotropies considered, 7 = 0.1 and 
7 — 0.5, the decay of QD with distance can be well 
fitted by an exponential function a + 6exp(— crt), 
where a, 5, and c are constants. Nevertheless, we 
notice that, while for A < 1 QD vanish exponen- 
tially, in the cases where A > 1, we obtain a con- 
stant long-distance value for QD that depends only 
on 7 and A. 

4. XXZ chain in the presence domain walls 

In order to consider QPTs of first-order and 
infinite-order (see, e.g., Ref. [13]), let us analyze 
the XXZ spin- 1/2 chain in the presence of a bound- 
ary magnetic field generating domain walls, whose 
Hamiltonian reads 4^, 49j 



JV-l 



(13) 
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Figure 2: Decay of quantum discord with the dis- 
tance between the spins sites for some values of A 
and 7. The points are the computed values of QD 
and the lines are the exponential fits (see the text 
for details). 
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Figure 3; (Color Online.) Behavior of QD as a func- 
tion of the distance for the XXZ chain with do- 
main walls for TV = 22 sites and magnetic field 
h/J — 5.0. The fitting curves in the gapless re- 
gion — 1 < A/J < 1 display the form Q{n) — 
a + bn^'^ (power-law decay) whereas in the an- 
tiferromagnetic and kink-type regions the form is 
Q{n) = a + b exp(— cn) (exponential decay), with 
a, b, and c constants. 



where the coupling J and A are exchange parame- 
ters, with the effective magnetic field h > repre- 
senting the interactions of the spins with the bound- 
ary surfaces. In order to focus the discussion on the 
values of the anisotropy and the magnetic field, let 
us assume, without loss of generality, J = 1. For 
A < — 1, the model is in an antiferromagnetic phase 
while for —1 < A < 1 the model is in a disordered 
critical (gapless) region. The model exhibits an 
infinite-order QPT at the antiferromagnetic point 
A = — 1 and a first-order QPT at ferromagnetic 
point A = 1. Moreover, for A > 1, as shown in 
Ref. [i^, the model presents a further first-order 
QPT governed by the critical field 



(14) 



which separates two quantum phases for a chain 
of arbitrary length: a ferromagnetic ground state 
{h < he) and kink- type ground state (h > he). 
Remarkably, he provides the exact location of the 
(first-order) phase transition for chains of any size, 
the critical field remaining fixed as the number of 
sites is changed. 

For the Hamiltonian given by Eq. (IT^ . we nu- 
merically compute QD for a finite-size chain by 
exact diagonalization. In the regime of tempera- 



ture T = 0, we take the system in its ground state. 
From the ground state wave-function, we then ob- 
tain the two-spin reduced density matrices for ar- 
bitrarily distant pairs. Since the system exhibits 
Z2-symmetry, the numerical evaluation of Eq. ([3]) 
can be performed similarly as in the XY model. 
Results for a chain with N — 22 sites is shown in 
Fig. [21 where QD is taken for the two-spin density 
operator pon, with pon standing for the pair {N/2, 
N/2 + n). Remarkably, a power-law decay is ob- 
served in the gapless region — 1 < A < 1 while 
exponential decay occurs for the antiferromagnetic 
and kink- type phases (for the ferromagnetic phase, 
QD is vanishing 36|, |37[). This observation keeps 
unchanged for distinct magnetic fields h as long as 
h keeps the system in a fixed quantum phase. 

5. Conclusion 

In summary, we have investigated the decay of 
pairwise quantum correlations (as given by quan- 
tum discord) as a function of distance along a spin 
chain. We have found distinct behaviors for the de- 
cay of QD as a function of distance as QCPs are 
crossed in the quantum phase diagram. In par- 
ticular, the onset of macroscopic order has been 
shown to be accompanied by the development of 
long-distance QD in the XY model. For the XXZ 
chain with domain walls, we have shown that the 
critical region exhibits power-law decay for QD 
whereas the gapful (antiferromagnetic and kink- 
type) phases display exponential decay. With re- 
gard to the fact that the set of zero-QD states has 
zero volume in the state space [Uj , we observe that 
the possible existence of quantum discord for dis- 
tant sites does not imply a priori in a power-law de- 
cay or in a decay to a constant value of QD. In fact, 
this is an interesting and surprising result, which is 
in contrast with the behavior of pa irwise entan- 
glement, which is typically short-ranged (see, e.g., 
the discussion in Ref. [j] and the Figs. 2(e) and 



2(f) in Ref. |35|). This characterization of QPTs 



in a condensed matter system from the point of 
view of information theory is an important aspect 
of the cross-fertilization of these two fields, turn- 
ing out to be possible only as we consider measures 
of quantumness that go beyond the entanglement- 
separability paradigm. In particular, a promising 
result presented here is the asymptotically constant 
behavior of QD in the XY model. Indeed, since QD 
has been recognized as a resource in several con- 
texts, it is potentially relevant to further investigate 
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the possible role of this long-range behavior of QD 
in quantum communication protocols [5(5] and its 
possible effects in biological systems i51|]. Further 
appealing directions for future research are the in- 
vestigation of the critical aspects of QD (correlation 
length, critical expo nents, etc.) and the behavior 
of multipartite measures of quantum cor relations 
(see, e.g., Ref. [EilHilHI) at QPTs. 
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Appendix A. Magnetization density and 
correlation functions 

The transverse magnetization density is given by 

4il 



where 



(1 + A cos 0) tanh {(3uj^) 



004, = (7Asin0)^ + (1 + Acos( 



d0, (A.l) 



(A.2) 



and (3 — l/ksT with ks being the Boltzmann's 
constant and T the absolute temperature. 



The two-point correlation functions read [42 1 



and 
where 



Go G_i 

G„_2 G„_3 

Gi Go 
G2 Gi 

Gn Gn-l 



G-n 
G-n+1 

G_i 

G-,i+2 
G-n+3 



(A.3) 



Gi 



— GnG- 



(A.4) 



(A.5) 



, , tanh(/3a;0) f a\ n ^ \ 
dcp — — {cos [nq)) (1 + A cos ( 



2TTUJd 



-7A sin (n0) sin (f)} 



(A.6) 
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